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PACS numbers: ????????? 02.50.-r, 05.40.-a, 05.70 Ln, 82.20-w 



key words: current fluctuations - Bethe ansatz - Symmetric exclusion 



I. INTRODUCTION 



Understanding how currents of particles or of heat fluctuate through non equilibrium systems has motivated lots of efforts over 
the last two decades. A number of general symmetries of the distribution of these fluctuations, based on the microreversibility 
of the dynamics, have been discovered ifTll |20T. 2l|_ 29||. Beyond these symmetries, the distribution of the current fluctuations 



has been calculated in several cases [22 J25 J26 , 27, 3§]. For diffusive systems, a general theory, the macroscopic fluctuation 
theory, has been developed 0] which, under some conditions 0,0, 53], a U° ws one t° calculate the whole distribution of these 
current fluctuations for systems maintained in a non-equilibrium steady state by contact with two reservoirs (of particles at 
different chemical potentials or of heat at different temperatures) 53 H3l- F° r driven systems, which are not diffusive, the 
fluctuations belong to the universality class of the Kardar Parisi Zh ang equ ation ll3"lll and can be related to the fluctuations in 
growth models and to the theory of random matrices lf30L F35L [36l [40t |42L |48L 

For diffusive systems as well as for driven systems, it has been very helpful to analyse simple models, in particular exclusion 
processes f32l l34ll47ll . In exclusion processes, one considers particles hopping on a lattice with a hard core interaction which 
prevents two particles from being on the same lattice site. In one dimension, exclusion processes can be solved by the Bethe 
ansatz |]4jll2lll3jllS|23|,[24, 39, 48l |49j |50J and a number of exact results have been obtained on the fluctuation of the current or 
the distribution of particles for a given initial condition lfT9l |j7i [38l l44ll . 

In the present work we consider the symmetric simple exclusion (SSEP) on an infinite one dimensional lattice. By definition 
of the model each lattice site is either empty or occupied by a single particle. The dynamics is stochastic: each particle hops to 
each of its neighboring sites at rate 1 if the move is not forbidden by the exclusion rule (each site is occupied by at most one 
particle). At time t = each site at the left of the origin (x < 0) is occupied with a probability p a and each site at the right of the 
origin (x > 0) is occupied with a probability pt, (see figure 1). We also assume that the measure at t = is Bernoulli, meaning 
that there are no correlation between the occupation numbers of the different sites in the initial condition. We call Q t the total 
flux of particles between site and site 1 during the time interval t. Our goal in the present paper is to determine, for large times 
t, the probability distribution of Q t . 
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FIG. 1: The initial condition when p a = 1 and pb — 

The same initial condition has already been considered for the asymmetric simple exclusion (ASEP) MM 5251 Si 
l49l[50ll . in particular in connection to random matrix theory. 
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Our main result is that, for large t, the generating function of this total flux is given by 

( e >«*) ~ eVSFC") (1) 

where the function F(u>) is defined by 



1 (\n+l i roc 

F(lu) = -=V u n = ~ dfclog 



1 + we~ fc " 



(2) 



and uj is a function of p a , p;,, A 

w = Pa(e A - 1) + p b (e- A - 1) + PaPb (e x - l)(e- A - 1) . (3) 



Remark: It has already been noticed 111411 that in the steady state of the SSEP on a large but finite chain, the generating function 
of the current was also a function of p a , pb, A through the same single parametrer, ui. This is a general property of the symmetric 
exclusion process that we shall establish in section 2 for a general graph. 

Remark: From (I1I2I31 > one can obtain all the cumulants of Qt, in the large t limit. They all grow like \ft and the 
prefactor in the cumulant (Q™) c is a polynomial of degree n in p a and pb- For example, for p a = p and pb = 0, one gets 



lim t-^(Q 2 t ) 

t— >oo 

lim t-V 2 (Q*) 
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while for p a = Pb = P one gets for the even cumulants (all the odd ones vanish by symmetry) 

lim t-^ a (QH) c = -*p(l-p), 

lunt-^(Qt) c = p(l-p)-^^(l-p)». 

t^oo -./7T ^/7T 



This can be compared to the variance of the position X t of a tagged particle Q 111 [3^14^1431 111]], which also grows like y/i, by 
arguing that the typical distance between particles is p^ 1 so that (Qf) ~ p 2 (Xf). 

From the knowledge of the generating function (Q]), one can obtain the distribution of Q t . This distribution takes, in the long 
time limit, a scaling form 

Pro ( % = ~ exp[VtG(g)] , (4) 



Vt 

where G{q) can be related to F(u>) by a Lagrange transform. This allows one to obtain the asymptotics of G(q). One can see 
from the integral representation © of F((J) that, for large positive ui, 

F(w) ~ ^ [log^] 3 / 2 + | [logo;]- 1 / 2 + .. 

and this gives for large positive q 

G(q)~-^q 3 + q\og(p a (l-p b )) + .. (5) 

This paper is organized as follows: in section II we prove that the generating function ^e A< 2*) of Q t is a function of the 
densities p a , Pb and of A through the single parameter u> defined in (f3j). In section III we obtain an exact expression of the 
generating function of Qt from which we derive the large t asymptotics (1112b . 
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FIG. 2: The large deviation function G(q) and its asymptotics for large q (dashed lines) in the cases p a — pb = 1/2 (left curves) and p a = 1, 
pb — (right curves). 



II. EVOLUTION OF THE CORRELATION FUNCTIONS AND THE w DEPENDENCE 



Let t x be a binary variable which indicates whether site x is occupied (t x = 1) or empty (t x — 0). In this section we are 
going to argue that all the ri-point correlation functions of these t x 's are polynomials of degree n in p a and This will allow 
us to show that the n-th moment (Q™) is also a polynomial of degree n in p a and pt- Then using the particle-hole symmetry, we 
will show that the generating function (e x Q l \ depends on p a , p^ and A through the single parameter u defined in (01. 

One can write the exact evolution of the expectation (t x ) where (.) denotes an average over the stochastic evolution and the 
initial conditions. On the infinite line it takes the form 

^(t x ) = (t x+1 ) + (t^x) - 2(t x ) . (6) 

As at time t = 0, one has (t x ) = p a for x < and (t x ) — pi, for x > 0, it is clear that at any later time the solution of © is 
linear in p a and pi,. 

Similarly for any x < y the evolution of the 2-point correlation function is given for y > x + 1 by 

^( T * T y) = ( T ^+l T !/) + { T x-\Ty) + (T x T y+ l) + (T x T y -l) - i{T x T y ) (7) 



and for y = x + 1 by 



J t (r x r x+1 ) = (t x _ iTx+1 ) + (r x r x+2 ) 2(r x r x+1 ) . (8) 



These evolution equations do not involve higher correlation functions and as at t = the only possible values of (T x T y ) are 
Pa, p a pb or p\, one can easily see that at any later time (r x T y ) remains a quadratic function of p a and pb- 

One can generalize this property of the SSEP to all equal-time or unequal-time correlation functions on a general graph. Let 
us consider the symmetric exclusion process on such a graph : the vertices are either empty or occupied by a single particle and 
each particle can hop at rate one to every empty site to which it is directly connected by the graph. Therefore if X = {x±, ..x n } 
is a set of n different sites on the graph and if (jx) denotes their 7i-point correlation function 

(Tx) = (r xl -T Xn ) , 

it is easy to see that the evolution of (tx) is given by 

%^ = E X>(y>**)Kw\{x»}> - fa)] (9) 



4 



where x( x > v) — 1 or indicates whether the edge (x, y) is present or not on the graph. If initially all the r x 's are uncorrelated, 
some sites being occupied with probability p a and all the other sites occupied with probability pb, all the equal-time rt-point 
correlation functions are obviously polynomials of degree exactly n in p a and pb at time t = 0. This property is preserved 
by the evolution (O. A similar reasoning allows one to show that all the unequal time n-point correlation functions are also 
polynomials in p a and p& of degree at most n. 

On this general graph let us call A the subset of sites occupied at time t = with probability p a (all the other sites being 
occupied with probability pb). Obviously one has (r^) = p a if x € A and (t x ) — pb if x ^ A. If Q t is the total flux out of the 
subset A during time t, one has 

Qt = XVx(o)-T S (t). (io) 

It is clear from ( TTOb that (Q") can be expressed in terms of unequal-time n-point correlation functions and is therefore a 
polynomial of degree n in p a and pb (if one or more sites are repeated more than once in the correlation function, this may give 
a polynomial of lower degree as t x = t x for all k > 1). 



To show that (e x< ^* ) depends only on the reduced parameter u), let us consider a finite lattice of N sites, where initially all the 
Na sites of the subset A are occupied with density p a , and all the remaining TV — Na sites are occupied with density p\,. One 
can write the generating function ^e A ^*) as 

N A N-N a 

(e AQt ) = E E P P a(l-P«) NA - p pUl-p b ) N - NA - q e^ x R p , q (e x ), (11) 

p=0 q=0 

where R Ptq (z) is a polynomial of degree p + q in z which depends on time and on the choice of the subset A. The reason for 
the factor e~ qX R p . q {e x ) is that this term corresponds to situations where initially there are p particles in the subgraph A and q 
particles in the rest of the graph so that Q t can only take the values Q t — —q, —q + 1, ..p. If one expands this expression in 
powers of p a and pb, one gets 

N A N-N A 

<e AQt )=E E tiple-^S^), (12) 

p=0 q=0 

where S Pl q(z) is also a polynomial of degree p + q. For (Q") to be a polynomial of degree n in p a and pb one needs that for A 
small 

S p , q (e x ) = O(A^) , 

so that the polynomial S PiQ has the form 

S P ,q{e x ) — s p ,g(e A — l) p+q , 

where s v<q is a number which depends on time, on the graph and on the subgraph A but does not depend on A. This implies that 
(fTTT i can be rewritten as 

N A N—N A 

<e AQt ) = E E s ™ " ^ t^( e " A - = G M eA ~ V'Pbk'* - !)) ■ < 13 > 

p=0 rj=0 

This shows clearly that ^e^') is already a function of only two reduced variables: p a {e x — 1) and Pb(e~ x — 1). 

Let us now use the particle-hole symmetry. In the SSEP, holes have exactly the same dynamics as particles. Therefore the 
generating function of the flux Q t is left invariant by the particle-hole symmetry (p a , pb,X) — > (1 — p , 1 — pb, — A). In terms of 
the function G of two variables defined in (13[ , this means that, for any p a , pb, A 

G(p a (e x - l), Pb (e- x - 1)) = G((l - p Q )( e - A - 1), (1 - p b )(e x - 1)) , (14) 

and so, if a = p a (e X — 1) and (3 = Pb{e~ x — 1) , one has 

G(a, 13) = G{e~ x - 1 + ae~\ e x - 1 + (3e x ) . (15) 
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As ( TT3T > is valid for any A, one can choose e A = 1 + (3, which leads to 

G(a, 0) = G(a + P + a/3,0) . (16) 

This completes the proof that G(a, (3) and therefore (e A< ^*) are functions of the single variable 

a + (3 + afi = p a (e x - 1) + p b {er x - 1) + p a (e x - l)ph(e~ x -1)=oj. 

Remark: This ui dependence of (e A< 5*) was already noticed for the SSEP on a finite lattice of length L connected at its 
extremities to two reservoirs of particles in ifTill . where it was, however, only obtained in the large t and L limit. A consequence 
of the above discussion is that the u> dependence remains valid for any system size L, at any time t, in more complicated 
geometries (in particular in higher dimensions), and for more complicated connections to the two reservoirs at densities p a and 
p b . To illustrate this claim, let us consider as in ifTill the SSEP on a one dimensional lattice of L sites, with particles injected on 
site 1 and L at rates a and S and removed from these two sites at rates 7 and (3. It is well known that these rates correspond 
to site 1 being connected to a reservoir at density p a = and site L to a reservoir at density p b = j^g- Instead of these 
input rates, one could think of site 1 being connected to a large number Ni of sites (which mimic the left reservoir), which 
are all initially at density p a , with an exchange rate (a + j)/Ni between site 1 and each of these Ni sites of the reservoir. 
Similarly one can replace the rates (3 and S by a large reservoir of Nl sites, initially at density p b = g^g, with an exchange 
rate ((3 + S)/Nl with site L. If inially all sites i < io are at density p a and all sites i > io + 1 are at density p b , we are in a 
situation where all the sites of the graph composed by the one-dimensional lattice and the reservoirs are initially occupied with 
probability either p a or p b . In this geometry the flux Q t is then simply the total flux of particles between site io and site io + 1, 
and therefore, for this flux, we know from (TBI that <^e A( ^ f ) depends on p a , p b and A through the single parameter uj. 

Because the generating function (e^Q* ) depends only on the single parameter to defined in (0, it can be written as 

( eAQt )=E S n(t)" n - (17) 

n>0 

Therefore, to determine the coefficients s„(i), one can limit the discussion to the particular case p b = where the analysis is 
easier. Then 

<e^) = E s »W^( eA - 1 )"- < 18 > 

n>0 

This can be viewed as an expansion of ^e A< ^*) in powers to p a . We are now going to express the coefficients s n (t) in terms of 
properties of exclusion process on an arbitrary graph. For simplicity we consider a finite graph. Initially, only the sites of the 
subset A are occupied, with probability p a : very much like in ( fTTT i one can write the generating function (e A< 3*) as 

(e^>= E^'(l-Pa)^- |B| £>o t (£ )(7 )e^, (19) 

ECA q=0 

where Prot(£', q) is the probability that q particles have escaped from A during t, given that at t = only its subset E was 
occupied by particles. Comparing ( fT~8T > and ( fT9] l, we see that the coefficients s n (t) can be expressed in terms of the probabilities 
Prot(-E, q). In fact, this relation takes a simple form, which can be easily understood by taking simultaneously the limits p a — > 
and A — > 00, at fixed p a e x , in (TT~8T > and ( fT9] l : this leads to 

s n (t)= P™ t (E,n), (20) 

ECA,\E\=n 

where the sum is over all the subsets E of n sites and Pro t (_E, n) is the probability that, if initially all the n sites of E are 
occupied, the rest of the graph being empty, all the n particles have escaped from A at time t. 

III. THE BETHE ANSATZ 

The evolution equations of the correlation functions, on the infinite line, can be solved via the Bethe ansatz : for a fixed initial 
configuration, i.e. {t x (0)} = {i] x }, it allows one to obtain the expression of all the correlations functions at any later time. For 
example the solution of (O is 

(Ty(t)) =J2 P t {1 \y\x)Vx, (2D 

X 
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where 



P t (1 \y\x) 



l*i|= 



dzi 
2iirzi 



y—x 

z{ exp 







(zi - 


-1-2V 




*1 / . 



(22) 



where the integration contour is a counterclockwise circle of radius r (for convenience in what follows we choose r <C 1). 
Similarly one can show that the solution of J7|8b is, for y\ < 2/2, 

(2), 



(T Vl {t)Ty 2 (t)) = ^ P t {Vl,V2\xi,X 2 )r]x 1 Vx 2 , 



(23) 



where 



Pt 2 \yi,y2\xi,x 2 ) = f 

J\zi 



dz\ dz2 



ziz 2 + 1 - 2z 2 
z l z 2 + 1 — 



x exp 



1 1 

z ± -\ h z 2 H 4 }t 

z\ z 2 



(24) 



In fact (fJTJ and (|23l can be extended to arbitrary correlation functions to give 



(T yi (t)Ty 2 (t)..Ty n (t)) = ^ Pt n) (yiiV2, ■■,yn\xi,X 2 , --Xn) Vx 1 Vx 2 --Vx n 

Xl<X2<--<X n 

where P t ' (yi, y 2 , ■■, Z/nl^i, £2, ■ ■%n) has a Bethe ansatz expression which generalizes d22l241 > (see AppendixlAland j48ll ). 
One can then use d20b to establish that the coefficients s n (t) are given by 

Sn(t) = ^ P t (VUV2, -,yn\Xl,X 2 , ..X n ) . 

xi<X2<-.<x n <0<y 1 <..<y n 



(25) 



(26) 



For example one gets from ( 1221 l24b for n = 1 and n = 2 



| Zl |=r 2l7T (1 - Zlf 



exp 



zi + 2 )t 

z\ 



(27) 



and 



*a(t) 



dz\ dzi 

\ Zl \ = \z 2 \=r 2l7TZl 2i7TZ 2 

z\zi 



.2 2 



(l-zi)(l-z 2 )(l-ziz 2 ) 2 
z x z 2 + 1 - 2z 2 s 



x exp 



(1 - Zi) 2 (l - ZiZ 2 ) 2 V^lZ2 + 1 - 2zi 

By symmetrizing the above expression over z\ and z 2 one can show that this can be rewritten as 

1 



1 1 

Z\ H hz 2 H 4 i 

2l ^2 



S2 (t) = i/ ^l^det 
2 /i Zi i =r 2z7r 2i7r 



ZfeZ; +l — 2zij 1 



exp 



L<fc,Z<2 



1 1 

21 + hz 2 H 4 t 

z\ z 2 



(28) 



(29) 



This expression can in fact be extended to arbitrary n, using an approach which follows closely recent works by Tracy and 
Widom [48, 49]. The derivation, which is detailed in AppendixlAl yields 



(e A *)=5>« Sn (*)=£ 



\z k \=r 



n 

,fe=i 



2i7T 



det 



ZfcZ; + 1 - 2zi J 1 



(30) 



L<fc,2<r 



It is known l45ll (see also eq.(7) of J49ll ) that for a Fredholm operator K which transforms a function / into a function 

Kfhy 

Kf(z)= I ^K(z,z')f(z>), 

J\z\=r 2 «7T 
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one has 



det(J + uoK) = 



„> J\zk\=r [ fc=1 



n 



2l7T 



det (K(z k ,zi)) 1<k l<n . 



This implies that 

log[det(J + ujK)} = tr[log(7 + wif )] = - ^ 
Comparing ( f30b and OTb . we see that by choosing 



n>0 



\zk\=r 



n dz k 



Lfc=l 



, , x exp (z+±-2)i 

K(z, z r ) = — ^ 2 Li. 

V ' y zz' + l-2z 



one gets d32l that 



log(e 



AQt 



E 

ra>l 



n 

,fc=i 



2iTT 



n 

fc=i 



exp 



2 i 



ZfcZfc+i + 1 - 2z fc ^ 



(with the convention that z n+ i = z\ in the n-th term). 

In Appendix B, we derive the large t behavior of the integrals in the r.h.s. of (|34| |. which leads to 

log(e^) ~ - V t^l [± = VtF( w ) , 



(31) 



(32) 



(33) 



(34) 



completing the derivation of (1112b . 



Remark : from the expressions (|2~T1 l22l 123124-1 l25l IAU . one can in principle calculate arbitrary correlation functions of 
the occupation numbers Tj at time t. For the one-point and the connected two-point functions one gets when p a — 1 and pb = 



(r x ) 



dze t(z+l/z-2) z x 



and (tvt„ 



dzdz f e t{z + l/z+z' + l/z'-4) z x Z ly 



2iirz l-z " J 4ir 2 zz' zz' + 1 - 2z' ' 

Then, carrying out an asymptotic analysis similar to the one performed in Appendix 151 yields for n — 1 and n = 2 



with 



Gx{X) = — 

V 7T 



X/2 



e~ u du and G 2 {X,Y) 



2V2tt 



e"" du . 



(35) 



(36) 



The scaling form ( 135136b of the two-point function is very reminiscent of what is known for the steady state of the SSEP 
on an open interval [46] of length L , where the two point function scales like L _1 and the n-point (connected) correlation 
function scales like L 1_ ™ 1 18]. It would be interesting to know whether the scaling form ( f35l > remains valid for n > 2. If so one 
could then try to determine the scaling functions G n in order to obtain the large deviation function of an arbitrary density profile 



IV. CONCLUSION 



In the present work, we have shown that, for a step initial density profile, the generating function of the integrated current Q t 
is a function of a single parameter u>, defined in (|3}, which takes a simple closed expression (Q]i. This ui dependence is also valid 
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for an arbitrary graph. In one dimension, for this non-equilibrium initial condition, the distribution of Q t is clearly not Gaussian 
with a tail which decays faster than a Gaussian ©■ 

It would be interesting to generalize our results to other diffusive systems and to try to calculate the distribution of Q t with 
the same step initial condition, in order to see under which conditions one could recover the non-Gaussian decay (0. The most 
promising approach, at the moment, is to try to generalize the macroscopic theory of Bertini, De Sole, Gabrielli, Jona-Lasinio, 
and Landim JaSSl t0 a non-steady state initial condition. 

A possible extension of the present work would be to introduce a weak asymmetry (of order t~ x l 2 ) in the hopping rates to 
understand the cross-over in the current fluctuations between the SSEP and the ASEP. 

Another interesting question would be to determine the scaling form (f35T > of all the higher correlation functions and to obtain 
the large deviation function of an arbitrary density profile, when the initial condition is, as in figure Q] far from a steady state 
situation. 

APPENDIX A: DETERMINANT EXPRESSION (30) OF (e XQt ) 

In this appendix we derive the general expression (f30b of the (s„) from (|26T >. Our derivation relies heavily on results obtained 
by Tracy and Widom on the particle trajectories of an ASEP model with the same geometry ll48l l49ll . In the present appendix we 
will use the Bethe ansatz itself ( B48I1 . Theorem 2.1), 



Pt n \yi,V2, ..,y n \x 1 ,x 2 , ..x n ) = ^sgn(cr)^" 



n 

,fc=i 



dzk 

2%'KZk 



t(z k +l/z k -2) Vk-Xa(k) 
6 °(k) 



n 

k<l 



z k zi + 1 - 2z k 



(Al) 



as well as two algebraic identities ([48], eq. (1.6) and I49B , eq. (7)) 



E s s n ( C7 ) 



,fc=l 



k<l 



n -. 



n 



Zk - zi 



Zkzi + l- 2zk 



z\-z n 



= det 



n*(i-**) 



z k zi + l-2z k J 1<kl<r 



(A2) 



(A3) 



In order to obtain s n from d26l i. one has to sum Pj (j/i, y2, .., ynl^ii x 2i ■•%n) over x i < ■ ■ < x n < and < y\ < .. < y v 
: this yields 



E IK 

xi<..<x n <0 k 



i 



Z\..Z n 



Lfc=i 



z\..z k 

- Z\..Zk 



and £ n 



v(fe) - n i 

0<y!<..<y n k fc=l 



Za(k)--Z a (n) 
- Z a (k)--Za(n) 



Therefore 



E 



sgn(cr) 



n 



dz k e t{zk + 1/zk 2) Z x ..Zk Z a (k)~Z<r{n) 



2mz\ 



- fc 1 - Zi..Zk 1 - 2o-(fe)..-2o-( n ) 

We are now in position to use iA2i , which yields 



n 



J Lk<l 



Zg(k)Zcr(l) + 1 - 2z CT ( fc ) 
Z k Zi + 1 - 2z k 



-q dz k e t ^ +1 / Zk -^ z x ..z k 



Lfc=l 



2mz k {l - z k ) \ — Z\..Zk 



n 

k<l 



Zl - Z k 



z k zi + l-2z k 



In order to eliminate the remaining x z _} z '*\ - , it is necessary to use (IA21 > again. As the factors in dA2b are of the form y£ 
we first relabel Zk — > z n _ k , so that 



" dzke^+^-V z k ..z n 

2i7TZfe(l - Zfc) l-Zfe..Z n 



n 



- z k 



n dzke t( Zk + l/ Zk -2) 

2l7TZ fe (l - Zfc) 



J Lfc>* 

ru<i(**-*i 



Z/cZz + 1 - 2zfc 



JJ(z fc z ; + 1 - 2z fc ) 
.fc<; 



rii 

k=l 
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We then replace the expression above by its average over all permutations of the (zk) in order to apply ( IA2t . The first term of 
the product is unchanged and the second term gets a sign factor sgn(cr). Hence 



E 

(7 

1 
;?! 



sgn(cr) 



« dzke t(z k + l/z k -2) 

2i-KZk{l - Zk) 



Y[( z k - zi) 

k<l 



T\k<l( z <?(k)Za(l) + 1 - 2z CT (fc)) 

UkM ZkZi + 1 ~ 2zfe ) 



k=l 



z a(k)-- z a(n) 
— z a{k)-- z a{n) 



Lfc=l 

Applying (IA3b then leads to 



which is the expression ( f3Qb . 



» dzfce t(^ + l/ Z ,-2) 

H 2i7r(l - zk) 2 



n 



ZfeZ; + 1 - 2zk 



Su(t) = —. 



ndzk 



Lfc=l 



dct 



e t{z k + l/z k -2) 

z k zi + l-2z k J 1<kl<n 



as 



APPENDIX B: DERIVATION OF THE ASYMPTOTICS @4} OF THE GENERATING FUNCTION 

In this appendix we derive the large t behavior of the integrals /„ = tr K n , which appear in J32I33I341 ) . I n can be expressed 

/„ = tvK? = <t> — —K t (z 1> z 2 )---K t (z n ,z 1 ) = <b [[ — 



2tTT 2z7T 



fe=l 



2iir z k z k +i + 1 - 2z k 



where the integration contour of each z k is a circle of radius r <§; 1 (and by convention we have z n+ k = Zfe). By using the 
identity 



t(z k + 1 + l/z k -2) i 
= |- / rfifc e t fc (z fc + l + l/z fc -2) ^ 

+ l/z fe - 2 Zfc+i + l/zfc-2 7 



/„ can be expressed as 



Sc{l,..,n} Lfc=l 



n 



2iirz k 



n 



1 



z fe+ i + l/z fe - 2 



n / ^ e 



tfc(zfc + l + l/z fc -2) 



If E ^ {1, .., n}, then there exists at last one k such that k ^ E. For fixed values of the (zi)i^ k and the (ti)i^ k , the integral 
over Zfc can be written either as 



r,tk-lZ k 



dz k 



2inz k z k +i + l/z k - 2 



or as 



1 



dz k 

2mz k (z k +i + l/z k - 2) (z k + l/zk-i - 2) ' 



As the integration contour is a circle of very small radius these expressions obviously vanish. Therefore the case E = {1, .., n} 
gives the only non-zero contribution to /„ : 

ft ^e^^ 



\z k \=rJ0 k 



2iirz k 



The integrals over the z k 's can now be evaluated using the saddle-point method around z k = yik/tk^i- This yields 



JJ dtk P -(V*T-v^T) 2 



k=l 



2y/TTt k 



l/2 r l 



f T\da k e- t{ak - ak - l) 
Jo 



Finally, since the integrand above is only non-vanishing when \a k — ai\ ~ for t — > oo, one gets 

%/2 r i 



J n ~ - / dai da 2 ---da n e- t ^ ak - a >'-^ 2 = J— - 
V 71 / Jo J-oo V nn 
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and this completes the proof of ( [34-b . 
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